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Experimental results on heavy-quark production in proton- proton (p+p), deuteron-gold 
(d+Au), and gold-gold (Au+Au) collisions at ysjvjv = 200 GeV at the Relativistic Heavy Ion 
Collider (RHIC) are reviewed. 
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Hadrons carrying heavy quarks, i.e. charm or bottom, are important probes in high energy 
hadronic collisions. Heavy quark-antiquark pairs are mainly produced in initial hard scattering 
processes of partons. While some of the produced pairs form bound quarkonia, the vast majority 
hadronizes into particles carrying open heavy flavor. The latter do not only provide a crucial 
baseline for quarkonia measurements but are also of prime interest on their own. Heavy-flavor 
measurements in p+p collisions provide an important proving ground for quantum chromody- 
namics (QCD). Because of the large quark masses, charm and bottom production can be treated 
perturbatively (pQCD) even at small momenta^. This is in distinct contrast to the production 
of particles carrying light quarks only which can be evaluated within the pQCD framework only 
for sufficiently large momenta. Systematic studies in p+p and d+Au collisions should be sensi- 
tive to the nucleon parton distribution functions as well as nuclear modifications of these such 
as shadowing^. In Au+Au collisions, heavy quarks present a unique probe for the created hot 
and dense medium. Imp ortan t observables in addition to heavv flavor yields ar e energy loss^^ 
and azimuthal anisotropyl^El as we ll as quarkonia suppression^ or enhancement^^!!!]! 

At RHIC, the PHENIX and STAR experiments study heavy-quark production indirectly via 
the measurement of electrons from semileptonic decays of hadrons carrying charm or bottom. 
In addition, STAR directly reconstructs hadronic D meson decays in p+p and d+Au collisions. 

The mid rapidity electr on spectrum from heavy-flavor decays measured by PHENIX in p+p 
collisions at ^/s = 200 GevEl 

is shown in Fig. ^ (left panel) . A le adin g order PYTHIA calcula- 
tion and a recent next-to-leading order (FONLL) pQCD calculation are compared to the data. 
The corresponding STAR measurement^^ agrees with these data within the substantial uncer- 
tainties. While the spectrum is significantly harder than predicted by PYTHIA, the FONLL 
calculation describes the shape better but still leaves room for further heavy-flavor production 
beyond the included NLO processes, e.g. via jet fragmentation. Bottom decays are expected to 
be essentially irrelevant for the electron cross section at px < 3 GeV/c and become significant 
only for px > 4 GeV/c. 

D° meson spectra, measured by STAR in p+p and d+Au collisions at ^/snn = 200 GeV are 
shown in Fig. ^ (right panel) together with electron spectra from heavy- flavor decays It is 
important to note that the electron and D meson data are compatible with each other, i.e. the 
measured electron spectra agree within errors with spectra calculated for semileptonic D meson 
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Figure 1: Invariant differential cross section of electrons from heavy- flavor decays in p+p collisions at 200 GeV 
in comparision with PYTHIA and FONLL pQCD calculations (left panel). Transverse momentum spectra of D° 
mesons and electrons from heavy-flavor decays in p+p and d+Au collisions at 200 GeV (right panel). 



decays. The nuclear modification factor RdA is calculated as the ratio of the d+Au electron 
spectrum to the spectrum from p+p collisions scaled with the number of underlying nucleon- 
nucleon binary collisions. Averaged over the range 1 < pr < 4 GeV/c, STAR measure 
RdA = 1-3 ± 0.3 ± 0.3, which is consistent within errors with binary scaling as expected for 
a point-like hard pQCD process. A modest Cronin enhancement is not excluded by the data. 
Electron spectra measu red by PHENIX in d+Au collisions as function of centrality confirm the 
observed binary scaling^] 




Figure 2: Heavy-flavor electron pr spectra for different Au+Au centrality selections at 200 GeV (scaled for clarity) 
compared to fits to the binary scaled p+p measurement (left panel) . Nuclear modification factor Raa in the range 
2.5 < pr < 5 GeV/c as function of the number of participant nucleons for tt° and electrons from heavy-flavor 

decays (right panel). 



For Au+Au collisions at ^/snn = 200 GeV, PHENIX has demonstrated^! that the yield 
of electrons from heavy- flavor decays is consistent with binary scaling in the range 0.8 < pr < 
4 GeV/c, which is entirely dominated by charm decays. Preliminary electron data, however, 
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Figure 3: Elliptic flow strength «2 of electrons from heavy-flavor decays in Au+Au collisions at 200 GeV as 
function of pr in comparision with recombination model calculations with and without charm quark flow. 



indicate a strong modification of the spectral shape in central collisions. Relative to binary 
scaled p+p data, electrons from heavy-flavor decays ar e sig nificantly suppressed at high px in 
central Au+Au collisions as shown in Fig. (left panel)'^. This observation is consistent with 
a scenario where quarks suffer energy loss while propagating through the hot and dense medium 
created at RHIC. The nuclear modification factor i?^ -5 , defined as the ratio of the yield of 
electrons from heavy-flavor decays in Au+Au collisions in the range 2.5 < pr < 5 GeV/c to the 
binary scaled yield in p+p collisions in the same pt range, is shown as function of the number 
of participant nucleons N part in Fig. |2] (right panel) together with the same quantity for neutral 
pions as measured by PHENIXHSI. The high pt electron suppression is comparable to the pion 
suppression, but the experimental uncertainties are still substantial and do not allow to establish 
the centrality dependence of heavy quark energy loss. 

A complementary observable related to the interaction of heavy quarks with the medium 
created in Au+Au collisions is elliptic flow. The elliptic flow strength v i fo r electrons from 
heavy-flavor decay s is shown as function of pt in Fig. El The PHENIX dataP3 and preliminary 
results from STAR^I are compared with two recombination model calculations^. Charm quark 
flow is consistent with the experimental data but the uncertainties are currently too large to 
exclude a scenario where charmed hadrons acquire a non-zero V2 through being produced via 
the coalescence of a non-flowing charm quark and a flowing light quark. 

PHENIX has measured J ftp in the dielectron channel at mid rapidity and in the dim uon 
channel at forward and backward rapidities in p+p and d+Au collisions at 200 GeV^^". The 
shape of the rapidity distribution shown in Fig. 0] (left panel) for p+p collisions agrees well 
with Color Octet Model (COM) and PYTHIA calculations using different parton distribution 
functions. The measured total cross section is consistent with predictions from COM and Color 
Evaporation Model (CEM) calculations. In d+Au collisions weak cold nuclear matter effect are 
observed. The ratio of the rapidity distribution from d+Au collisions to the properly scaled p+p 
measurement show in Fig. |1] (right panel) indicates both weak absorption as well as weak shad- 
owing of the gluon distribution function in nuclear matterl^. Given the current uncertainties, 
it is difficult to disentangle these small effects. 

In a deconfined medium the yield oiJ/ip might either be reduced due to the expected screen- 
ing of the att ractiv e QCD potential^ or possibly even enhanced via coalescence^ or statistical 
recombinationEE3 "While scenarios lea din g to a strong J/ijj enhancement are disfavored by first 
data from Au+Au collisions at 200 a definite answer will emerge only from the currently 

ongoing analysis of the high statistics Au+Au data sample recorded in RHIC Run-4. 
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Figure 4: J/ip differential cross section, multiplied with the dilepton branching ratio, as function of rapidity in 
p+p collisions at 200 GeV compared with a PYTHIA calculation (left panel). Ratio of d+Au and appropriately 
scaled p+p J/V> rapidity distributions compared with model calculations 21,22 (right panel). 
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